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ABSTRACT 

Motivated  by  recent  observations  of  solar  hard  X-ray  bursts  with  a 
spatial  resolution  of  a  few  arc  seconds,  we  compute  the  theoretically  pre¬ 
dicted  spatial  variation  of  hard  X-ray  flux  versus  height  in  both  thick 
target  non-thermal,  and  thermal,  models  of  solar  hard  X-ray  bursts.  Our 
work  on  the  thick  target  model  revises  previous  results  in  this  area  by 
adopting  a  more  realistic  model  for  the  flaring  atmospheric  structure,  and 
also  by  taking  into  account  energy  loss  and  scattering  processes  in  the 
evolution  of  the  non-thermal  electron  beam  which  have  previously  been 
neglected. 

It  is  pointed  out  that  in  the  so-called  "thermal”  model  currently  in 
vogue,  there  is  a  substantial  non-thermal  bremsstrahlung  component  at  all 
photon  energies;  it  is  further  shown  that  this  non-thermal  component  results 
in  a  hard  X-ray  flux  versus  height  distribution  whose  maximum  brightness 
per  unit  length  is,  similar  to  the  thick  target  model,  located  in  the 
chromosphere,  and  not  in  the  thermal  source  itself  (where  the  bulk  of  the 
total  low  photon  energy  bremsstrahlung  is  emitted) .  The  characteristics  of 
the  hard  X-ray  height  distributions  in  both  thick  target  and  thermal  models 
are  summarized  and  compared,  and  predictions  as  to  the  observable  spatial 
structure  of  the  X-ray  flare  in  both  cases  are  made. 


Subject  Headings:  Sun:  chromosphere  -  Sun:  corona  -  Sun:  flares  -  Sun:  X-rays 


I.  INTRODUCTION 


Current  theoretical  models  of  the  hard  X-ray  bursts  commonly  associated 
with  solar  flares  fall  into  two  main  classes  -  thermal  and  non-thermal.  In 
the  former  model  the  X-rays  are  attributed  to  collisional  bremsstrahlung 
within  an  ensemble  of  thermally  relaxed  electrons  (e.g.  Chubb,  Kreplin,  and 
Friedman  1966;  Brown  1974;  Crannell  et  al.  1978;  Matzler  et  al.  1978;  Brown, 
Melrose  and  Spicer  1979;  Smith  and  Lilliequist  1979;  Emslie  and  Brown  1980; 
Emslie  and  Vlahos  1980;  Smith  and  Auer  1980;  Emslie  1981a).  In  the  latter 
the  X-rays  similarly  result  from  bremsstrahlung,  although  this  time  from  a 
beam  of  suprathermal  electrons  incident  on  a  cold  background  target.  The 
exact  target  geometry  involved  has  led  to  a  variety  of  non-thermal  models  — 
thick  target  (Brown  1971;  Emslie  1980),  thin  target  (Datlowe  and  Lin  1973), 
and  trap  (Takakura  and  Kai  1966;  Brown  and  Hoyng  1975;  Melrose  and  Brown 
1976;  Hudson  1978;  Emslie,  McCaig,  and  Brown  1979). 

Current  X-ray  spectral  measurements  (e.g.  Hoyng,  Brown  and  van  Beek 
1976)  are  unable  to  distinguish  between  the  various  models  outlined  above, 
because  of  the  instability  of  the  X-ray  spectrum  deconvolution  problem  (e.g. 
Craig  1979)  and  because  of  the  ability  of  a  thermal  distribution  with  a  non- 
uniform  temperature  to  mimic  a  "non-thermal"  (e.g.,  power-law)  spectral  form 
(Brown  1974).  Similarly,  available  polarization  measurements  (Tindo  et  al. 
1970,  1972a,  b;  Tindo,  Mandel'stam,  and  Shuryghin  1973;  Tindo,  Shuryghin, 
and  Steffen  1976)  are  not  reliable  enough  (Mandel'stam  1980)  to  discriminate 
amongst  the  predictions  of  various  models  (Brown  1972;  Langer  and  Petrosian 
1977;  Bai  and  Ramaty  1978;  Emslie  and  Brown  1980;  Leach  and  Petrosian  1981). 
It  has  therefore  been  suggested  (e.g.  Rust  and  Emslie  1979;  Emslie  and  Rust 
1980)  that  spatially  resolved  hard  X-ray  images  may  provide  vital  clues  to 
aid  in  the  discrimination  amongst  candidate  models. 
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Observationally,  the  only  "spatially  resolved"  hard  X-ray  data  published 
at  the  time  of  writing  is  the  fortuitous  observation  of  a  behind- the- limb 
flare  by  two  satellites  separated  by  a  suitable  distance  in  ecliptic  longi¬ 
tude  (Kane  et  al.  1979).  However,  with  the  successful  launch  of  the  NASA 
Solar  Maximum  Mission  satellite,  with  its  Hard  X-Ray  Imaging  Spectrometer  (van 
Beek  et  al.  1980)  on  board,  there  is  now  an  abundance  of  spatially  resolved 
(£8"  resolution)  data  available  for  study,  and  preliminary  results  have  been 
presented  (de  Jager  et  al.  1981;  Hoyng  e£  al .  1981). 

Theoretically,  the  only  published  study  of  the  height  distribution  of 
solar  hard  X-rays  in  the  context  of  a  particular  model  is  by  Brown  and  McCly- 
mont  (1975),  who  analyzed  the  thick  and  thin  target  non-thermal  models.  How¬ 
ever,  the  predicted  distribution  in  these  models  is  a  sensitive  function  of 
the  density  structure  of  the  model  atmosphere  adopted;  Brown  and  McClymont 
(1975),  in  making  use  of  a  quiet  Sun  model  (Vemazza,  Avrett,  and  Loeser 
1973),  in  fact  obtain  results  which  significantly  differ  from  those  obtained 
by  using  a  more  realistic  flare  model  atmosphere  (Machado  and  Linsky  1975; 
Machado  et  al.  1980);  we  show  this  below  in  §111. 

In  the  present  paper  we  therefore  revise  Brown  and  McClymont's  (1975) 
calculations,  using  a  more  realistic  model  atmosphere  as  a  target,  and  also 
taking  into  account  the  effects  of  collisional  scattering  and  reverse  current 
ohmic  energy  losses  on  the  non-thermal  electron  beam  (effects  neglected  by 
Brown  and  McClymont  in  their  analytic  calculations).  We  also  compute  the 
height  distribution  of  hard  X-rays  in  the  dissipative  thermal  model  (Brown, 
Melrose,  and  Spicer  1979;  Smith  and  Lilliequist  1979;  Smith  and  Auer  1980), 
giving  due  recognition  to  the  non-thermal  high  energy  tail  component  in  the 
electron  velocity  distribution  in  such  models  (Brown,  Craig,  and  Karpen  1980; 
Emslie  and  Vlahos  1980) .  The  hard  X-ray  height  distribution  in  coronal  trap 
models  (Takakura  and  Kai  1966;  Melrose  and  Brown  1976)  depends  on  the  magnetic 
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field  geometry  In  the  trap  and  on  the  electron  scattering  process  within  the 
trap  (Melrose  and  Brown  1976;  Wentzel  1976),  and  will  not  be  dealt  with  here. 

In  §11  we  describe  the  method  by  which  our  (numerical)  results  are  obtained. 

In  §111  we  compare  our  thick  target  results  with  those  of  Brown  and  McClymont 
(1975).  In  § IV  we  study  the  thermal  model  and  show  that,  contrary  to  some 
expectations  (Rust  and  Emslie  1979;  Emslie  and  Rust  1980),  the  hard  X-rays 
in  this  model  do  not  have  their  brightest  component  in  the  thermal  part  of 
the  flaring  loop,  although  the  bulk  of  the  total  photon  flux  at  low  energies 
emanates  from  the  (larger)  thermal  region  .In  §V  we  discuss  the  results  obtained 
and  show  how  the  two  models  may  be  distinguished  by  spatial  measurements  with 
sufficiently  high  resolution;  such  spatial  resolution  is  within  current  instru¬ 
mental  capabilities. 
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II.  METHOD  OF  COMPUTATION 

Since  we  shall  assume  the  hard  X-ray  emission  from  any  thermal  source 
to  be  spatially  uniform  (see  §IV),  it  remains  only  to  describe  the  computa¬ 
tion  of  the  non-thermal  component  present  in  both  the  thick  target  and 
dissipative  thermal  models.  We  shall  make  use  of  the  procedure  developed 
by  Emslie  (1978,  1980)  to  compute  the  dynamics  of,  and  X-ray  bremsstrahlung 
yield  from,  the  beam  of  non-thermal  electrons;  we  note  that  his  mean  scattering 
and  energy  loss  theory  yields  results  which  are  very  similar  to  a  full  Fokker- 
Planck  analysis  of  the  beam-target  interaction  (Leach  and  Petrosian  1981). 

It  has  been  pointed  out  by  Chambe  and  Henoux  (1979)  that  Emslie' s 
(1978)  treatment  of  ion-neutral  collisions  neglects  the  inelastic  momentum 
loss  from  the  charged  beam  particles.  This  loss  is  given  by 


Inelastic 


1  d£ 

mV|j  dt  ’ 


(1) 


where  \||  is  the  component  of  velocity  in  the  precollision  direction,  and 
m  and  E  the  beam  particle  mass  and  energy  respectively.  Adding  these  in¬ 
elastic  momentum  loss  terms  to  Emslie' s  (1978)  equation  (22)  gives  the 
revised  total  momentum  loss  rate  in  a  target  of  arbitrary  ionization  level: 


In  this  equation  m0  and  m^  are  the  electron  and  proton  masses  respective¬ 
ly,  e  (e.s.u.)  is  the  electronic  charge,  x  is  the  fractional  ionization 


number  of  hydrogen,  and  A,  A',  A"  are  Coulomb  logarithms,  defined  in  Emslie's 
(1978)  paper.  It  follows  that  Emslie’s  (1978)  equation  (25)  for  the  scatter¬ 
ing  parameters  8  should  read 


a  „  2*A  +  ( 1— x)  A”  .  « 
Pe  A’  +  x  (A-A’)  »  Pp 


(3) 


(Be  and  8^  corresponding  to  electron  and  proton  bombardment  respectively). 
In  turn,  equation  (14)  of  Emslie  (1980),  which  includes  the  effect  of  both 
collisions  and  reverse  current  ohmic  energy  and  momentum  losses  from  the 
beam  electrons,  should  read 
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where  p  is  the  direction  cosine  (w.r.t.  the  vertical)  of  the  beam  electrons, 

n  and  n  are  the  plasma  density  (cm-3)  and  resistivity  (e.s.u.)  respect- 
* 

ively,  Eq  is  the  minimum  injection  energy  required  of  an  electron  to  pene- 

—2 

trate  to  a  column  depth  N  (cm  ),  and  Fqo,  and  6  define  the  injected 

electron  flux  spectrum: 

F  (E  )  (electrons  cm  %  *keV  l)  ■  (5-1)  F  E^  1  (E  +E  )  ^  .  (5) 

o  o  oo  oo  oo  o 

Note  that  Emslie's  (1980)  equation  (13)  remains  unaltered.1 
*We  have  re-performed  Emslie's  (1980)  calculations  using  the  revised  version 
of  his  equations  (13)  and  (14)  and  find  no  significant  changes  from  his  results. 

Using  the  revised  equations  as  outlined  above,  we  may  now  derive  the 
variation  of  electron  flux  with  depth  in  the  target,  and  so  compute  the  X-ray 
bremsstrahlung  yield  as  a  function  of  depth  and  photon  energy  (see  Emslie's 
[1980]  equation  [22]).  Emslie  (1980)  integrated  this  over  depth  to  obtain 
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the  resulting  bremsstrahlung  spectrum  only;  however,  it  is  a  straightforward 
matter  to  use  his  results  directly  to  compute  the  height  variation  of  the 
X-ray  flux  at  a  given  photon  energy.  Results  follow  in  the  next  two  sections. 


.  i 
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III.  RESULTS  FOR  THE  THICK  TARGET  MODEL 
As  mentioned  in  §1,  the  earlier  treatment  of  this  problem  by  Brown  and 
McClymont  (1975)  uses  the  quiet  Sun  atmospheric  density  structure  (Vemazza, 
Avrett,  and  Loeser  1973)  as  a  target.  This  structure  in  fact  bears  little 
resemblance  to  that  found  in  flares  (Machado  and  Linsky  1975;  Lites  and  Cook 
1979;  Machado  et  al.  1980).  The  principal  difference  between  the  quiet  Sun 
and  flare  model  atmospheres  is  the  particle  column  depth  (cm  2)  of  the 
transition  zone,  with  its  associated  steep  density  gradient;  in  empirical 
flare  model  atmospheres  N  =1021cm  2,  while  N  is  -2  x  10iacm  2  in  the  quiet 

Let 

Sun  (Vernazza,  Avrett,  and  Loeser  1973)  and  -3xl019cm  2  in  active  regions  (Basri 
et  al.  1979).  It  is  not  immediately  obvious  whether  one  should  employ  an 
active  region  or  flare  empirical  model  as  a  target  for  the  bremsstrahlung- 
producing  non-thermal  electrons,  since  they  are  flash-phase  phenomena  and 
empirical  flare  modeling  to  date  (e.g.  Machado  and  Emslie  1979;  Machado  et  al. 
1980)  applies  only  to  the  later  stages  of  the  flare.  However,  since  in  a 
thick  target  model  the  principal  cause  of  the  increase  in  N  is  the  bombard- 
ment  by  the  bremsstrahlung-producing  electrons,  and  since  it  is  easily  shown 
that  preflare  chromospheric  material  rises  in  temperature  by  the  required 
107K  or  so  in  a  very  short  time  (see  eq.  [8]  of  Emslie  and  Noyes  1978) , 
it  seems  more  reasonable  to  use  the  flare  models  of  Machado  et  al.  (1980) 
as  a  background  target  for  our  computations.  (For  further  discussion  of 
the  validity  of  empirical  flare  model  atmospheres,  see  Emslie,  Brown,  and 
Machado  1981.)  Now,  the  predicted  height  distribution  of  thick  target  hard 
X-rays  is  a  combination  of  two  effects:  (i)  density  variation  (the  thick 
target  yield  is  proportional  to  the  ambient  number  density),  and  (ii) 
electron  flux  attenuation  (the  yield  is  also  proportional  to  the  electron 
flux  at  a  given  level,  which  is  in  turn  a  decreasing  function  of  the  column 
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density).  Thus  one  sees  that  by  moving  the  level  at  which  the  large  photon 
flux  enhancement  occurs  due  to  effect  (i)  to  a  larger  (flare)  value  of  the 
column  density  N  (and  so  to  a  level  at  which  the  electron  flux  has  attenuated 
to  a  much  higher  degree),  a  significant  difference  from  Brown  and  McClymont's 
(1975)  predicted  X-ray  flux  versus  column  density  profile  may  result. 

In  addition  to  these  considerations,  we  note  that  Brown  and  McClymont 
(1975)  neglect  collisional  scattering  of  the  beam  electrons,  and  also 
reverse  current  ohmic  energy  losses  from  the  beam;  these  processes  can 
result  in  substantial  modifications  to  the  electron  flux  versus  depth  profile 
(Emslie  1980).  Further,  Brown  and  McClymont's  (1975)  assumption  of  a  fully 
ionized  target  underestimates  the  bremsstrahlung  yield  at  greater  depths, 
where  the  fractional  ionization  level  of  hydrogen  is  small  (Brown  1973). 

All  these  effects  are  incorporated  in  the  present  treatment. 

We  have  calculated  the  variations  of  hard  X-ray  flux  versus  depth  for 
two  model  atmospheres.  The  first  uses  the  model  chromosphere  'FI'  of 
Machado  et  al.  (1980)  and  the  second  the  model  chromosphere  'F2'  of  these 
authors  (hereafter  these  models  will  be  referred  to  as  MAVN1  and  MAVN2 
respectively).  To  these  models  we  attached  isothermal  coronae  (temperature 
3  x  107  K)  3  x  104km  in  length  (measured  down  one  half  of  the  loop).  The  density 
at  the  loop  apex  was  taken  to  be  3  x  1010cm  3  for  model  MAVN1  and  10ncm  3 
for  model  MAVN2,  and  the  density  variation  along  the  coronal  part  of  the 
loop  was  taken  to  be  smoothly  increasing  with  depth  such  that  the  integrated 
column  density  N^z  at  the  transition  zone  matched  Machado  et  al. ' s  (1980) 
values  of  1.5xl02°  cm  2  and  1.6  x  1021  cm  2  respectively.  These  coronal 
models  are  admittedly  somewhat  crude;  however,  it  turns  out  that  the  coronal 
hard  X-ray  emission  per  unit  length  in  the  thick  target  model  is  negligible  compared 
to  that  in  the  chromosphere  (Figures  1  and  2),  so  that  any  deficiencies  in  the 
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the  model  adopted  will  not  significantly  affect  the  conclusions  of  the  paner. 

In  Figures  1  and  2  we  show  the  variation  of  <|>f  the  fraction  of  the  total 

X-ray  flux  emitted  per  km  of  vertical  distance  at  depth  z  (km)  (column  depth 

N  [cm  2]  ),  for  these  two  flare  models.  The  electrons  were  assumed  to  be 

injected  by  a  "black  box"  mechanism  (with  no  intrinsic  hard  X-ray  signature) 

at  the  top  of  the  coronal  loop  (N=0)  with  a  flux  spectrum  given  by  equation 

(5)  with  F  *  1019  electrons  cm  2s  1 ,  E  =20  keV,  and  <5  =  4  (see  Emslie 
oo  ’oo 

1980).  The  dashed  line  indicates  a  purely  collisional  calculation  of  the 
electron  beam  dynamics  and  the  solid  line  a  calculation  involving  reverse  current 
losses  as  well.  These  reverse  current  losses  are  quite  sensitive  to  the  atmospheric 
model  adopted,  particularly  in  the  corona  (Emslie  1980),  which  is  why  we  show 
the  purely  collisional  calculation  for  comparison,  even  though  it  is 
unphysical.  Note  also  that  the  requirement  that  the  reverse  current  be 
stable  to  current-driven  plasma  turbulence  sets  an  upper  limit  to  the 
significance  of  reverse  current  ohmic  dissipation  in  the  electron  beam 
energetics.  Only  for  certain  coronal  conditions  (e.g.  high  density,  approx¬ 
imately  equal  electron  and  ion  temperatures  )  is  such  a  large  flux  Fqo 

allowed  to  pass  stably;  if  these  conditions  are  not  met  then  F  must  be 

oo 

reduced  and  the  effects  of  reverse  currents  correspondingly  decreased 

(Emslie  1981b).  Note  finally  that  the  figures  only  extend  down  to  a  depth  N  - 

1022  cm  2.  This  is  the  thermalization  depth  of  electrons  with  initial  (i.e. 

injected)  energy  ^300  keV  (e.g.  Emslie  1978);  since  data  on  hard  X-rays 

with  photon  energies  of  this  order  are  poor,  we  consider  extrapolation  to 

higher  electron  energies  (and  so  greater  depths)  than  this  highly  uncertain. 

If  the  spectrum  (5)  does  in  fact  extend  to  arbitrarily  large  energies  E  , 

o 

then  4>(N)  will  asymptotically  approach  a  power-law  form  with  index  ( 1—6 ) /2 
=  -1.5  (Brown  and  McClymont  1975;  Emslie  1978). 
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Figure  2  It  is  evident  from  the  figures  that  the  effect  of  reverse  current 

ohmic  dissipation  is  significant  for  the  coronal  region  of  model  MAVN1 
only;  this  is  due  to  the  relatively  low  coronal  densities  and  consequently 
small  collisional  losses  (compared  to  ohmic  ones)  in  this  region  (Emslie 
1980).  Note  that  the  effect  of  the  reverse  current  is  to  enhance  the 
coronal  component  of  the  hard  X-ray  emission  relative  to  the  chromospheric 
component.  In  fact  (see  Emslie  1980) ,  both  chromospheric  and  coronal 
components  are  reduced  (in  absolute  flux  units)  from  their  purely  collisional 
values;  however  the  effect  is  more  severe  on  the  chromospheric  emission  due 
to  the  reduction  in  flux  penetrating  to  large  depths  effected  by  the  reverse 
current  ohmic  dissipation  in  the  corona. 

Comparing  our  results  to  those  of  Brown  and  McClymont  (1975)  we  notice 
that  in  model  MAVN2  (Figure  2)  the  region  of  brightest  X-ray  emission  per 
unit  length  is  narrower  than  in  either  Brown  and  McClymont 's  (1975)  results 
or  the  results  for  model  MAVN1  (Figure  1)  and  is  located  at  the  base  of 
the  transition  region,  as  opposed  to  mid-way  into  the  chromosphere.  This 
is  because,  as  remarked  above,  the  transition  zone  level  in  model  MAVN2 
is  so  deep  that  the  electron  beam  has  already  significantly  attenuated 
before  it  reaches  the  chromosphere,  to  such  an  extent  that  the  increase  of 
chromospheric  density  with  depth  from  the  base  of  the  transition  zone  down¬ 
wards  cannot  overcome  this  diminishing  flux  effect.  In  the  quiet  Sun  and 
MAVN1  models,  however,  this  coronal  attenuation  is  not  so  significant,  and 
consequently  the  region  of  maximum  bremsstrahlung  yield  per  unit  length  is 
pushed  downwards  into  the  chromosphere  and  made  broader. 

These  differences  are,  however,  probably  quite  undistinguishable  by 
observation,  at  the  current  level  of  instrument  capability.  The  height 
range  over  which  <f>  ^  ^’^max  is  ~2x103  kin  (3"  in  angle)  for  Brown  and 


McClymont's  quiet  Sun  calculation  (see  their  Figure  4),  =  103  kra  (14")  for 

model  MAVN1  (Figure  1),  and  =5xl02  km  (1")  for  model  MAVN2  (Figure  2).  All 
of  these  angular  sizes  are  below  the  currently  available  spatial  resolution 
of  hard  X-ray  detectors  (van  Beek  ££  al.  1980). 

It  is  worth  commenting  further  on  the  results  for  model  MAVN2.  Figure  2 
shows  that  the  depth  dependence  of  hard  X-ray  emission  in  the  coronal  part 
of  the  loop  is  significantly  different  for  low  (e  »  10  keV)  and  high  (e  -  70 
keV)  photon  energies.  For  low  energies  the  principal  cause  of  the 
decrease  of  4>  with  depth  is  effect  (ii)  above, namely  electron  flux  attenuation, 
while  for  high  energies  effect  (i)  (increasing  density  with  depth)  dominates, 
causing  the  X-ray  flux  per  unit  length  to  Increase  with  depth.  One  would  expect 
effect  (ii)  to  begin  to  dominate  effect  (i)  (causing  a  maximum  in  :)  when  one 
reaches  the  depth  at  which  electrons  of  energy  E  ’■>  e  begin  to  lose  energy  rapidl 
(note  the  dependence  of  the  Bethe-Heitler  bremsstrahlung  cross-section  [ Heitler 
1954J  on  electron  energy  E) .  This  occurs  when  (e  [kevj  )2  £  6»  e'AN  3 

(N/1017)  (Emslie  1978);  the  curve  for  e  -  40  keV  in  Figure  2  has  a  maximum 
at  N  -  3  x  1020  cm  2,  in  good  agreement  with  this  prediction. 

This  concludes  our  discussion  of  the  thick  target  model.  In  the  next 
section  we  shall  consider  the  thermal  model,  and  compare  the  results  from 


the  two  models  in  §  V. 
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IV.  RESULTS  FOR  THE  THERMAL  MODEL 

In  this  section  we  investigate  the  height  distribution  of  hard  X-rays 
expected  in  the  "dissipative  thermal"  model,  originally  proposed  by  Brown, 
Melrose,  and  Spicer  (1979),  and  currently  under  extensive  investigation  (see 
Smith  and  Lilliequist  1979;  Emslie  and  Brown  1980;  Emslie  and  Vlahos  1980; 

Smith  and  Auer  1980;  Smith  and  Brown  1980;  Emslie  1981a).  This  model  is 
characterized  by  a  "thermal"  distribution  of  electrons  at  the  top  of  a  coronal 
loop,  confined  parallel  to  the  loop  axis  by  a  strong  magnetic  field  and  per¬ 
pendicular  to  the  loop  axis  by  a  pair  of  collisionless  ion-acoustic  turbulent 
fronts,  which  result  from  the  instability  of  the  reverse  current  associated 
with  the  outward  heat  flux  from  the  source  (see,  e.g.  Kahler  1975)  and 
which  move  along  the  arch  at  approximately  the  local  ion-sound  speed  (Smith 
and  Lilliequist  1979;  Smith  and  Auer  1980). 

It  was  recognized  by  Brown,  Melrose,  and  Spicer  (1979)  and  by  Vlahos 

and  Papadopoulos  (1979)  that  in  such  a  model  very  high  energy  electrons  are 

not  confined  by  the  turbulent  plasma  fronts  and  escape  downwards  in  much  the 

same  manner  as  a  thick  target  beam.  Although  the  exact  details  of  this  escape 

process  are  somewhat  uncertain  (Smith  and  Brown  1980),  it  is  typically  found 

that  all  electrons  with  velocities  v  £  3  v  (where  v  is  the  mean  electron 

e  e 

thermal  velocity)  can  escape.  Brown,  Craig,  and  Karpen  (1980)  and  Emslie  and 
Vlahos  (1980)  have  considered  the  thick  target  bremsstrahlung  yield  from 
these  high  energy  electrons  and  find  that  it  is  comparable  to  the  thermal  yield 
for  photon  energies  e  greater  than  a  few  times  ^n^v2  (the  mean  electron 
thermal  energy).  Specifically,  Emslie  and  Vlahos  (1980)  give  the  ratio  of 
thick  target  to  thermal  bremsstrahlung  at  photon  energy  e  as 


1(e)  -  l 


6.6x10*  1  B4 

nL  4(5-2) 


e-B2/2  ge/kT  (kT)Z  ln 


(2glkT); 


e  <  E 


min 


(6) 


;  e  >  E 


min 


a 
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where  k  is  Boltzmann's  constant,  S  the  ratio  of  the  critical  escape 

velocity  to  vg  and  6  the  spectral  index  of  the  power-law  escaping  tail 

(see  Emslie  and  Vlahos  1980); note  that  assuming  the  tail  to  be  a  simple 

extended  Maxwellian  [see  Brown,  Craig,  and  Karpen  1980]  only  reduces  n  by 

a  factor  of  ^,2  for  the  parameter  set  considered  and  so  does  not  affect 

the  conclusions  to  follow  .  n  (cm’3),  T(K),  and  L(km),  are  the  density,  electron 

temperature,  and  length  of  the  thermal  source  and  is  the  threshold  energy 

for  an  electron  to  escape  through  the  conduction  fronts  ■  *j82kT.  We  shall 

hereafter  consider  a  source  model  with  8=3,  5  *  4,  n  *  1011  cm  3  (see  Smith 

and  Lilliequist  1979;  §  III)  and  T  =  22.2  keV  (so  that  the  electron  escape 

energy  threshold  E  is  100  keV) ;  results  for  other  values  of  these  parameters 
min 

.  i 

can  readily  be  obtained  using  equations  (6)  and  (8) .  For  the  parameter  set 
chosen, 

n(e)  »  - exp  (0.045  e)ln(^p.)  (7) 


Figure  3 


(with  e  measured  in  keV) ;  this  is  shown  graphically  in  Figure  3. 

If  we  assume  that  the  thermal  bremsstrahlung  emission  is  spatially  uniform 


over  z  <  z*  (N  <  N*) ,  then  the  fractional  (relative  to  the  combined  thermal  plus 
thick  target  non-thermal  emission  in  one  leg  of  the  arch)  yield  per  km  of  loop  is 


<$>(£»  z) 


2  +  n(e)<J>tc  (£,z)L 

(1  +  n(e) ]L 


n Ce)4>cc  (e,z) 
l  +  n(e) 


;  z  <  z’ 


;  z  >  z 


(8) 


where  $  (£,z)  is  the  thick  target  fractional  yield  as  found  in  §111,  but 

with  a  lower  cutoff  of  100  keV  introduced  to  the  spectrum  (5). 

Figure  4  jn  Figure  4  we  show  the  resulting  <Kz)  structures  for  various  e;  results 

for  the. thick  target  model  neglecting  reverse  currents  (see  Figure  2)  are  shown 


for  comparison.  Model  MAVN2  was  used  in  Che  calculations;  the  results  using 


model  MAVN1  are  similar.  The  column  density  N*  was  taken  to  be  lO20  cm 
and  the  corresponding  source  length  2X10  cm  (note  the  symmetry  in  the 
source  model) .  Note  that  the  ratio  of  thick  target  to  thermal  values  of 
$  in  the  model  is  in  fact  independent  of  L  :  to  see  this  we  rewrite  equa¬ 
tion  (7)  as 


-2 


n(e) 


n  (e) 


(9) 


with 


n  <e)  -  369  exp  (0.045  e)ln  — )  * 

°  €  ! 


(10) 


Then  equation  (8)  becomes 


#(e,  z) 


2  +  no(e)<j>tt  (e,z) 
L  +  nQ(e) 


;  z  <  z* 


n0<e)4>tt  (e,z) 

L  +  nQ(e) 


;  Z  >  z* 


(11) 


so  that  the  ratio  of  $  values  in  the  'Chromospheric” and  "thermal"  parts  of 
the  source  is 


9(e) 


V£)t*tt 

2  +  Vemtt 


^£ * ^chromosphere 
^G,Z^thermal  region 


(12) 


Since  for  all  values  of  e  the  second  term  in  the  denominator  is  negligible 
compared  to  unity  (the  suprathermal  electrons  give  off  little  bremsstrahlung 
in  the  relatively  tenuous  coronal  plasma  -  note  the  sharp  drop  off  in  t  at 
z>  z*  -  Figure  4),  we  see  that  0(e)  is,  to  a  good  approximation, constant : 
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9(e) 


Vemtt  ^£* chromosphere 
2 


(13) 


It  is  immediately  evident  from  Figure  4  that  the  maximum  bremsstrahlung 
yield  per  unit  length  at  all  energies  still  occurs  at  the  base  of  the  transi¬ 
tion  region.  This  is  contrary  to  earlier  speculations  (e.g.  Rust  and  Emslie 
1979;  Emslie  and  Rust  1980)  that  the  most  Intense  hard  X-ray  emission  should 
be  emitted  in  the  thermal  region  at  the  top  of  the  arch.  Comparison  of  the 
thick  target  and  thermal  model  results  in  Figure  4  shows  that  the  fractional 
contribution  from  the  coronal  region  is  indeed  enhanced  in  the  thermal  model; 
however  this  is  not  at  the  expense  of  the  chromospheric  emission  but  rather 
at  the  expense  of  the  emission  in  the  region  between  the  conduction  fronts 
and  the  transition  region  (i.e.  between  z  ■  -2  x  104  km  and  z  *  0)i 
'Note  that  D  (z)  is  an  increasing  function  of  z  for  all  energies  in  this 
region  of  the  atmosphere;  only  when  N  reaches  -1021cm  2,  the  thermalization 

depth  for  electrons  of  energy  E  ■  100  keV,  does  any  significant  decrease  in  $ 

mm 

occur;  cf.  discussion  of  thick  target  model  curves  in  §  III. 

This  behavior  is  in  fact  easily  understood:  since  the  electron  energy  required 
to  penetrate  to  chromospheric  depths  in  model  MAVN2  is  80  keV  (see,  e.g., 
Emslie  1978),  the  fractional  chromospheric  yield  in  the  thick  target  model  is 
matched  in  the  thermal  model  by  the  contribution  from  the  high  energy  escaping 
tail.  Further,  due  to  the  high  background  plasma  density  in  the  chromosphere, 
the  high  energy  tail  thermalizes  in  a  very  short  distance  and  also  has  a  high 
bremsstrahlung  yield  per  unit  length  (the  additional  fact  that  the  chromo¬ 
spheric  plasma  is  only  partially  ionized  further  enhances  the  chromospheric 
emission  -  see  Brown  1973).  However,  in  the  tenuous  coronal  plasma  between 
the  confining  turbulent  conduction  fronts  and  the  transition  region,  the 
thick  target  component,  being  inefficient  compared  to  thermal  emission  at 
these  densities  (Smith  and  Lilllequist  1979;  Smith  and  Auer  1980)  is  very 


small.  We  thus  see  that  the  height  distribution  of  hard  X-rays  in  the  dissi¬ 
pative  thermal  model  is  characterized  by  tvo  very  low  intensity  regions,  one 
in  each  leg  of  the  flaring  loop,  flanked  above  and  below  by  relatively  bright 
regions,  making  five  regions  in  all.  The  more  intense  of  these  bright  regions 
are  the  lower  ones,  with  the  degree  of  contrast  between  them  and  the  uppermost 
region  increasing  with  increasing  photon  energy  (see  Figure  4) . 
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V.  DISCUSSION 

The  principal  result  of  the  above  analysis  is  that  the  brightest  hard 
X-ray  region  is  always  at  the  base  of  the  transition  zone,  in  both  thick 
target  and  thermal  models.  This  result  follows  in  the  thick  target  model 
because  of  the  use  of  a  more  realistic  model  atmosphere  than  previously 
employed  (cf.  the  results  of  Brown  and  McClymont  1975),  and  in  the  thermal 
model  because  of  the  presence  of  a  precipitating  high  energy  non-thermal  tail 
of  electrons,  which  behave  exactly  as  a  thick  target  beam. 

In  our  discussion  of  the  thermal  model  (§IV)  we  have  tacitly  assumed 
that  the  high  energy  tail  is  always  present  (i.e.  we  have  assumed  a  steady- 
state  situation) .  Processes  by  which  this  escaping  tail  may  be  repopulated 
from  the  bulk  thermal  plasma  are  not  well  understood  (seeEmslie  and  Vlahos  1980; 
Smith  and  Brown  1980)  and  so  it  must  be  conceded  that  perhaps  the  X-ray  height 
distribution  depicted  in  Figure  4  strictly  applies  only  to  an  initial  tran¬ 
sient  state.  In  such  a  case  the  "relaxed"  $  distribution  would  simply  be  a 
rectangle  of  length  L/2  and  height  2/L  (see  Figure  4).  Clearly  whether  or  not  the 
high  energy  tail  is  continuously  repopulated  is  of  crucial  importance;  for 
the  present  we  note  that  the  interpretation  of  microwave  bursts  by  Emslie 
and  Vlahos  (1980)  in  the  context  of  the  dissipative  thermal  model  demands 
that  the  tail  does  exist. 

Returning  to  the  results  of  Figure  4,  we  note  that  the  distribution  will 


in  fact  be  time-dependent,  as  the  conduction  fronts  advance  along  the  arch 

(recall,  however,  the  results  of  5  IV  [  equation  (12)  J, 

which  showed  that  this  variation  of  L  does  not  affect  the  dominance  of 


the  chromospheric  component  in  the  $ [z ]  profile) .  An  obvious  corollary 
to  this  is  that  the  source  has  only  a  finite  lifetime  T  %  L^/c^  ,  where 


Hz 


is  the  distance  from  the  energy  release  point  to  the  transition  zone 
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Jj  —1 

and  cg  is  the  ion  sound  speed  -  104  r  cn  s  .  For  plausible  parameters 
we  see  that  t  is  of  order  30  seconds,  which  is  somewhat  greater  than  the 
time  resolution  of  available  hard  X-ray  imaging  instruments  (van  Beek  et  al. 
1980).  Thus  it  appears  that  the  introduction  of  a  "multiple  kernel  model", 
such  as  has  been  advocated  by  Brown,  Craig,  and  Karpen  (1980)  to  explain 
observed  (Miitzler  e_t  al.  1978)  emission  measure/temperature  correlations  (see, 
however,  Emalie  1981a)  may  be  unnecessary.  (Note,  however,  that  even  if  such  a 
multiple  kernel  model  were  appropriate, with  the  kernels  all  formed  near  the 
top  of  the  arch,  we  should  replace  L  by  a  mean  L  in  the  analysis  of  §IV, 
thereby  resulting  in  a  distribution  very  similar  to  that  of  Figure  4  [recall 
che  insensitivity  of  the  results  to  L  —  § IV ] ) . 

Finally  we  note  that  the  best  spatial  resolution  currently  available 
in  the  hard  X-ray  spectral  region  is  rather  limited  (^  8"  ;  van  Beek  jst  2d. 

1980).  An  angular  resolution  of  this  order  corresponds  to  a  length  of 

*  5000  km  on  the  solar  surface,  implying  that  we  should  really  scan  Figures 

1,  2,  and  4  with  a  rather  broad  "filter"  of  this  width.  The  effect  of  this 
smoothing  depends  on  whether  we  are  considering  thermal  or  non-thermal 
models,  and  on  the  model  atmosphere  adopted.  For  instance,  for  the  thick 
target  model  with  atmosphere  MAVN1  (Figure  1),  one  coronal  "pixel"  will 
contain  some  10  5  x  5  x  103  =  52  of  the  emission,  while  a  chromospheric 
pixel  will  contain  some  70£  of  the  total  emission  (from  one  leg  of  the  arch). 
Thus  In  this  case  we  should  expect  to  see  two  chromospheric  "bright  points" 
at  the  feet  of  the  flare  loop,  with  little  emission  elsewhere. 

In  the  thick  target  model  with  atmosphere  MAVN2,  we  find  that,  due  to  the 
greater  amount  of  coronal  material  compared  to  atmosphere  MAVN1,  a  coronal 
pixel  now  observes  some  15%  of  the  emission  (at  c  ®  40  keV) ,  while  the 

chromospheric  pixel  observes  a  similar  fraction  of  the  total  emission.  Thus 

we  should  here  expect  to  see  almost  uniform  emission  over  the  whole  loop. 
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Finally,  in  the  thermal  model,  we  see  that  a  pixel  observing  the  thermal 
region  will  contain  some  30-40?  of  the  emission,  a  chromospheric  pixel  will 
observe  some  20-30%,  and  a  pixel  in  the  region  between  the  thermal  region  and 
the  transition  zone  will  observe  at  most  2  or  3  percent  of  the  emission.  In 
this  case  we  therefore  expect  to  see  the  characteristic  "five  layer"  structure 
mentioned  in  §  IV. 

These  observational  predictions  are  sufficiently  distinct  to  distinguish 
between  the  above  three  scenarios  with  currently  available  instrumentation 
(van  Beek  at  al.  1980) ,  and  it  is  hoped  that  results  from  the  Solar  Maximum 
Mission  Hard  X-Ray  Imaging  Spectrometer  will  be  able  to  cast  some  light  upon 
which  scenario  actually  applies  to  solar  flares. 
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FIGURE  CAPTIONS 

Figure  1:  Fraction  of  total  thick  target  bremsstrahlung  yield  t>  (km  *) 

versus  depth  z  (km)  and  column  depth  N  (cm  2)  for  photon  energies 
10,  40,  and  70  keV  with  model  atmosphere  MAVN1  (Machado  et  al_.  1980) 
as  background.  The  z -coordinate  is  measured  from  the  transition 
zone  downwards;  note  the  different  scale  in  the  coronal  and 
chromospheric  regions  of  the  figure.  The  solid  curves  include 
the  effects  of  both  collisional  and  reverse  current  ohmic  losses 
from  the  bremsstrahlung-producing  non-thermal  electron  beam; 
the  dashed  line  takes  into  account  only  collisional  losses. 

Compare  with  Figure  2  and  with  Figure  3  of  Brown  and  McClymont 
(1975). 

Figure  2:  As  for  Figure  1,  except  using  atmospheric  model  MAVN2  (Machado 

et  al.  1980)  as  background.  Since  MAVN2  has  a  deeper  transition 
zone  than  MAVN1  the  electron  flux  has  already  been  significantly 
attenuated  by  the  time  it  reaches  the  transition  zone;  thus  the 
region  of  maximum  bremsstrahlung  yield  per  unit  length  is  very 
localized  and  occurs  at  the  base  of  the  transition  zone.  This 
result  is  markedly  different  to  the  results  of  Brown  and  McClymont 
(1975),  due  to  their  adoption  of  an  unrealistic  (quiet  Sun) 
background  atmosphere. 

Figure  3:  Ratio  of  thick  target  to  thermal  hard  X-ray  bremsstrahlung  (n) 
as  a  function  of  photon  energy  £  in  the  "thermal"  model  dis¬ 
cussed  in  §IV.  Above  the  electron  escape  energy  threshold  of 
100  keV  there  is  insignificant  thermal  emission  and  n  becomes 
infinite.  Note  the  substantial  contribution  of  the  thick  target 
emission  at  all  energies;  this  has  important  implications  for 
the  height  structure  of  hard  X-rays  in  this  model  (see  §IV  and 
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Figure  4) . 

Figure  4:  <J>  versus  N  for  thermal  model  with  parameter  set  outlined  in  §IV. 

The  figure  is  normalized  so  that  yj  d  z  is  unity  for  one  leg  of 
the  arch.  The  hard  X-ray  fluxversus  height  structure  is 
characterized  by  two  dark  regions  (one  in  each  leg  of  the  arch) , 
flanked  outside  and  in  between  by  bright  regions.  Note  that  the 
<p  (N)  structure  is  time  dependent  due  to  the  motion  of  the  conduc¬ 
tion  fronts  through  the  corona,  although  this  does  not  affect  the 
relative  sizes  of  <1>  values  in  the  corona  and  chromosphere  (see 
§IV). 
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